Derivation of quantum master equation with counting fields by monitoring a probe 



Takaaki MonnaQ 
^Department of Applied Physics, Osaka City University, 
3-3-138 Sugimoto, Sumiyoshi-ku, Osaka 558-8585, Japan 
(Dated: October 27, 2010) 

We show a microscopic derivation of a quantum master equation with counting terms which 
describes the electron statistics. A localized spin behaves as a probe whose precession angle monitors 
the net electron current by the magnetic-moment interaction. The probe Hamiltonian is proportional 
to the current, and is determined self-consistently for a model of a quantum dot. Then it turns out 
' that the quantum master equation for the spin-precession contains the counting terms. As an 

s ! , application, we show the fluctuation theorem for the electron current. 
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I. INTRODUCTION 

Recently the method of generating function has been extended to count the particle number flowing in quantum 

S junctions based on the quantum master equation with counting fields [l|, Indeed, the scheme has been successfully 
applied to the electron transport in quantum dots[l|,y], which provides a kinetic description of the recent experiments 
such as bidirectional electron transport in double quantum dots[3|, Hjj- Then the electron counting statistics illustrates 
the stochastic trajectory 0, Q. However, it remains unclear how the counting terms are derived microscopically from 
the interaction between the current and probe, i.e. we should give physical justifications of counting fields in the 
quantum master equation. The problem here is to derive the quantum master equation with counting terms from the 
total Hamiltonian including the probe. Although we only concern with the dynamics of the probe, the elimination of 
the system degrees of freedom is not self-evident due to the quantum correlation between the system and probe [H, 0]. 
Also, inclusion of the probe amounts to a self consisitent determination of the total Hamiltonian as shown in Sec. II. 

In this article, we address this issue based on a spin probe for the electron transport [jj, and give a physical 
justification of counting terms which appear in the quantum master equation. As an application, we show that the 
O • electron current satisfies a universal symmetry of fluctuation theorem [JU S HHUI EH, UM- 

This paper is organized as follows. In Sec. II, we present a model of electron transport and an idealistic spin probe 
. which monitors the current. In Sec. Ill, a quantum master equation with counting terms is derived. In Sec. IV, we 
derive the fluctuation theorem for net electron transfer based on the quantum master equation. 
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II. MODEL 
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Let us consider the electron current between the left and right reservoirs [li 151 [Iol - [l2j . As a model of a quantum 
dot, we consider a subsystem with a discrete energy level located between the reservoirs In order to detect the 
magnetic field caused by the electron current, we prepare a spin probe sufficiently near the current Q. The total 
Hamiltonian including the probe is given as 



H = H + V~ ^a z I (1) 
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where 



Hq = Mla + a + J dkhukja^jCikj 



V = dkH(ukjaa^j + v kj a k ja + ) 

|[/ dka+ iakll H\. (2) 



j'=l 

7 = 



i?o consists of Hamiltonians of the subsystem and reservoirs, where a, a k \, and a k2 are annihilation operators of the 
subsystem, left and right reservoirs, respectively. In the case of electron transport, they satisfy the anti-commutation 
relations {a, a + } — 1, {a ki , a k 'j} — 5ij8{k — k'). Also V presents a bilinear coupling between the system and reservoirs, 
which yields the current. According to the Ampere's law, the current yields magnetic field. Then the spin-probe shows 
a precession whose precession angle is proportional to the current. The amplitude of the magnetic field is proportional 
to that of the current /, and the magnetic moment interaction energy is written as — 4per z J. Here A is a measure 
of the coupling strength between the current and spin-probe, and a z is z component of Pauli matrices. The particle 
current I is defined as the time derivative of the particle number of the left reservoir. It is also possible to concern 
with the particle number of the right reservoir 

I' =^[J dka+ 2 a k2 ,H], (3) 

which is equal to I in the stationary state. To concern only with the left reservoir means that the probe monitors the 
electron moving between the left reservoir and the subsystem. 

Note that the definition of the total Hamiltonian H includes the current J, which depends on H. This condition 
amounts to a self-consistent equation for H, which is solved in the following subsection. 

A. Self-consistent determination of the total Hamiltonian 

In this subsection, we construct the total Hamiltonian which is consistent with the definition of the current I defined 
in Eq.Q. As the lowest order evaluation of the current, let us define Iq as 

h = ™[ / dka+ iakl ,Ho + V). (4) 

It is straightforward to verify that 



1 
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Then the first order evaluation of the total Hamiltonian is 



I = -— I dkh{u k \aa^ x - v kl a k ia + ). (5) 



H 1= H Q + V -^a z I . (6) 



Similarly, the next order evaluation of the current is given as 

h = ~Jfrlf dka^aki, Ht] 



2 " z {ihf 

And the second order expression of the total Hamiltonian is 



= h + -7T a z7-ky2 I dkh(u k iaa+ 1 +v k ia k ia + ). (7) 



H 2 = H Q + V -^a z h. (8) 
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In this way, we can recursively define the m-th evaluation of the total Hamiltonian H m as 

Xh 

H m — H() + V —<J z Im-l 

Im = dka kl a kl ! H m ] . (9) 

Then it turns out that linin^oo H n satisfies the self-consistent equation 

A/i 1 f 

lim H n = H + V + —<J z — [ dfca+Ofci, hm H n ]. (10) 

n— >oo Z 111 J n— >oo 

More explicitly, we have 

lim H n 

n— >oo 



H + j dkh{u k i{X)aa'l 1 + v kl {X)a k ia + ) + J dkh{u k2 aa\ 2 + v k2 a k2 a + ) 
1 + —a z 

y kiW = 7a — v ki- ( n ) 



B. Time evolution of the probe 



Depending on the value of a z , we use the abbreviated notation 

#A = (t| lim Hnlt); 

n— >oo 

H_A = UI lim H n \i). (12) 

n— >-oo 

Suppose that the initial state of the total system is described by the 2x2 density matrix 

p(0) = ( Pn Pn ) . (13) 

Then the unitary time evolution is described as 

Since we only concern with the precession angle of the probe, the trace is taken for the subsystem and reservoirs 
variables. Then the diagonal elements are invariant, while the off-diagonal elements evolves as 

Tr s , r e-*^* m (0)e* H -^. (15) 

As pointed out in Ref.@, Q, this quantity is identified as the characteristic function of the precession angle. In the 
following section, we pursue this issue in term of the corresponding quantum master equation. 



III. QUANTUM MASTER EQUATION 

In this section, we show our main result, i.e. a derivation of the quantum master equation with counting terms. 
Let us assume the weak coupling for V between the subsystem and reservoirs [3, 0]j an d pursue the dynamics of the 
off-diagonal element of the density matrix. As an initial condition, we assume that the left and right reservoirs are 
in mutually different equilibrium states and described by grand canonical ensembles, and the total density matrix is 
given as 



(16) 
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where p s (0) is the density matrix of the subsystem, and p r is that of the reservoirs 

Here /3y and /ij are the inverse temperature and chemical potential, and Sj is the grand partition function of the 
j-th reservoir. Then a quantum master equation for the system state Tr r /? is obtained in the Markovian limit: In the 
second order perturbation with respect to V, the relaxation time of the subsystem t s is evaluated as 

^ = <K, (18) 

where u is a measure of the coupling strength between the subsystem and reservoirs, and r c is the typical correlation 
time of the reservoirs. In the weak coupling limit, and for sufficiently short r c , we have 

T S > T C , (19) 

which justifies the Markovian description. 

Let us derive the quantum master equation. For this purpose, we use the interaction picture given as 

Pn(t) 

= e-i Hot a{t)e* Hot . (20) 

Then the von Neumann equation is written as 

Wf {t) = ^ HotV We-* Hot ,cT(t)] x , (21) 
where we have abbreviated the interaction Hamiltonian plus the probe Hamiltonian as 

V(X) =H X - H . (22) 

Also the generalized commutator is defined as 

[V(X),A] X = V{X)A-AV(-X) (23) 
for an arbitrary observable A. We solve the equation as 

<t) = ff(0) + 4 f ds ^ V ^ s )' a ( s )^> ( 24 ) 
in J 

where V(X,s) = ei H ° s V(X)e^i H ° s . By substituting Eq.flHJ into the right- hand- side of Eq.([2Tj). von Neumann 
equation is expressed as 

d , \ 

= -]hj ds l V ( X > t )dV(X 7S ) 7 a( S )] x } x 

+ hv(X,t),a(0)] x . (25) 
in 

Note that 

V(X,t) = J2[ dkh(u kj (X)aa+ j e t( - n ~^ t + v kj (X)a kj a + e- i( - n - u "^), (26) 

3 

where aside from Uki(X) and Vki(X) defined in Eq.(ll), we present Uk2(X) = Uk2 and ^2 (A) = Ufe2- Then for sufficiently 
long time, the Integral f* dse 1 ^-^ 1 " 8 ^ is replaced by / °° dse^-"^-'^ = ij^- + nS(Q - uj kj ), where P 
and S(x) denote the principal value and Dirac-delta. It is also remarked that in the weak coupling limit, a(t) in the 
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Eq. (l25l is replaced by Tr Y a(t) ® p r . By taking the trace over the reservoir variables, and calculating the expectation 
values with respect to the number state {|0), |1)}, one obtains the quantum master equation for Tr(n|cr(i)|n) as 



Tr r tT n (<) \_f -hi ho \ ( Tr r CTii(t) 
Tr r (T 00 (t) / I fcoi -fc o / I Tr r (TQo(t) J ' 



(27) 



where the transition rates are given as 



hi = Ai ^ X2 ri(/3i(m - Hi)) + A 2 r](l3 2 (hn - fi 2 )) 
Cl ° = Al ( T—i\ ) v(Pi(t& - Mi)) + A 2 i](l32(hn - fi 2 )) 




fcoi = A 1 | ) (1 - ^(Wl - mi))) + A a (l - T}(fo(m - &))) 

fcoo = Ai— ^(1 - - Mi))) +A 2 (1- V (f3 2 {h{l - M2 ))). (28) 

Here we abbreviated as A, = J" dfc27r|itfcj| 2 (>(il — uikj)- 

Since the interaction with the probe A should be sufficiently small to avoid the effect of measurements, the coefficients 
of quantum master equation (|28p is well-approximated as 

1 

" 1 



2 



1 ± — 

2 , 



= e^. (29) 



With this approximation by exponential e , it becomes now clear that the off-diagonal element of the density 
matrix (TT51) indeed counts the number of electrons, i.e. the coefficient of the phase factor e lXrn of Eqs. (l28l29p gives 
the probability P{t 1 m) that the net number of electrons moving from the left reservoir within time t is just m: 

2 

£ Tr r c^(i) = £ P(t, my Xm . (30) 

3 = 1 m 

Also, the present scheme shows a mechanism to include the counting terms e ±lA to the quantum master equation, 
which has been just assumedfj, d, H, [l2l [l3|. 



IV. FLUCTUATION THEOREM 



In order to further make clear that Eqs. (115127 1) indeed provide the electron counting statistics, we derive the universal 
symmetry of fluctuation theorem [lrH, l8l-ll(L_ Il3j . This symmetry contains several out of equilibrium relations 



such as linear and non linear response theory 10 , T3| . Since we have successfully derived the quantum master equation 
with counting terms, the fluctuation theorem is derived in a similar way as Ref. ;lj. A different point is the presence 
of the imaginary unit in the exponent of Eq. (|29l) . Let us calculate the eigenvalues of the matrix 



-fcn 
koi 



ho 
-fcoo 



(31) 



where we have assumed the first order approximation 



fc u = A x r){(3 x {m - mi)) + A 2V {(3 2 (hrt - M2 )) 

fcio = Axe-^rjipxiKl - mi)) + A 2V {p 2 {Kl - &)) 

fcoi - A ie tX {\ - ri(Px(m - mi))) + A 2 (l - n{fii(m - m 2 ))) 

fcoo =A l (l- 7](Pi(Ml - mi))) + A 2 (l - 77(/3 2 (^0 - m 2 ))). 



(32) 
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Then the eigenvalues of the matrix (jBlj) is evaluated as 
£f±(£A) 



1/ . . ^ , /Mi + 2 



2 (A 1 + ^ 2 )±^ 2 z j -A^aai-e^Xl-jjiJ^ + Cl-e-^Cl -»»)), (33) 

where we have abbreviated ry^ = rj(/3j(hQ — fJij)). The eigenvalues satisfy a symmetry 

g±{pi{m - mi) - 0a(Jifi - M2) - »A) = <7±(iA). (34) 

Let us derive the fluctuation theorem for the probability P(t, m) of the electron transfer in the long time regime. 
The Fourier series expansion of the probability P(t, m) in Eq. (|30[) is expressed by a linear combination of the normal 
modes 

J2 p {t, m)e lXm = ai e 9+(lA)t + a 2 e 9 - {lX)t , (35) 

m 

where the coefhcients are given as 

«1(2) 

1 1 ^ -(+)A 1 (2 m -l) + A 2 (2r l2 -l) \ 

V 4 V - 4^ ((! - elA )(! - m)V* + (! " e- u )r?i(l - »&))/ 

where the constants Ci^ are determined from the initial state. 

Without loss of generality, suppose that Reg + (iX) > Reg- (iX) . Then Eq. (|3"5j) amounts to 



lim ilogVP(t,m)e lAm = ff+ (zA) 



rn 

= g+{pi{m - mi) - /3 2 (M^ - M2 ) - iA) 

= lim ilogV p( t m ) e Wi(MJ-/*0-Ai(«i-p a )-<A)rr. (37) 
m 

where Eq. {[34]) is used in the second line. Thus, we have for the finite time 

J2P(t,m)e lXm = Ci(A,i)e 9 +( iA )* 

m 

J2p(t : m)e {f3l{rin '^ ) ' Mm ^^ ) ' iX)m = C 2 (X,t)e 9+{lX)t , (38) 

m 

where the amplitudes Ci(A, t) behave slower than exponential e ff +^ A '* as a function of time. By inverse transformations 
of Eq. (|3"8)) . we have the fluctuation theorem for the net electron transfer in the long time regime 

lim log ^f' m) . = m %(Hn - /ia) - - Ml)) . (39) 

which states that the ratio of the probabilities P(t, ±m) for the net number of electrons flowing through the probe 
should be balanced only by the thermodynamic affinity [lOl HT|. 



V. SUMMARY 



As a model of junction systems, we considered the electron transport between the subsystem, and two reservoirs. The 
electron current is monitored by the magnetic moment interaction with a localized spin probe. The total Hamiltonian 
is determined self consistently, since the total Hamiltonian includes current which depends on the Hamiltonian. Then 
we derive the quantum master equation with counting fields which describes the dynamics of the spin precession. 
This procedure gives a microscopic justification of the quantum master equation with counting fields, i.e. kinetic 
description of the electron transport. As an application, we have verified the fluctuation theorem for the electron 
transfer. 
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